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ABSTRACT 

We present the results of an optical multi-band (BVR) photometric monitoring program of 
22 core-dominated radio-loud quasars (CRLQs) and 22 radio-quiet quasars (RQQs). The aim 
was to compare the properties of microvariability in both types of quasars. We detected optical 
microvariability in 5 RQQs and 4 CRLQs. Our results confirm that microvariability in RQQs 
may be as frequent as in CRLQs. In addition we compare microvariability duty cycles in different 
bands. Finally, the implications for the origin of the microvariations are briefly discussed. 
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1. INTRODUCTION 

Optical Microvariability (OM hereafter) de- 
fined as flux changes on timescales ranging from 
minutes to hours, and with amplitudes that range 
from a few hundredths to a few tenths of a mag- 
nitude, was reporte d almost simultaneously as 
the quasar discovery (|Matthews fc Sandage lll963t 
Okelll967h . However, those reports were ignored 
or skeptically attributed to instrumental and/or 
weather causes. Nowadays, with the CCD tech- 
nology, and better statistical techniques, there is 
no doubt about the reality of OM. 

The importance of OM studies resides on 
the opportunity that they provide to set lim- 
its on the size of the emitting regions, to con- 
strain emission models, and to probe physical 
conditions in the innermost regions of active 
galactic nuclei (AGN), very close to the central 
black hole. Some OM studies have been car- 
ried out to investigate its incidence in radio quiet 
quasars ( RQQs) as compared to r a dio loud quasars 



(RLQs) dJang fc Miller 
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Stalin et al.ll2005l ; ICarini et al. l200"7llde Diego et al 



19981 . hereafter Paper I). The motivation for 
searching OM in RQQs and compare its inci- 
dence with OM in RLQs, is to identify the source 
of OM and try to constrain the most probable 
models: accretion d i sk an d /or jet emission (e.g. 
Lightman fc Earlev I Il974t iBlandford fc Konigl 



1991 



1979 1: iMarscher fc Gear I Il985l: iGeorge fc Fabian 



iQian et al 



1991; 



Marschcr ct al 
Gopa -Krishna et al. ll993a|bl:lMangalam fc Wiita 



19921: 
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19931 : IChakrabarti fc Wiita1l994l : lKrishan fc Wiita 
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Some of these studies have found that OM oc- 
curs in quasars regardles s of their radi o prop erties. 
The results reported by IStalin et all (|2004l ) indi- 
cate that although strong radio emission does not 
guarantee microvariability generation the most 
variable objects are those with mo re intense ra 
dio emission. On the other hand, ICarini et al. 



(|2007l ) found a higher instance of OM in radio- 
loud quasars than in radio quiet quasar. Accord- 
ing to these results, orientation could play a fun- 
damental role. Thus, objects with jets aligned 
closer to the line of sight are those with larger 
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duty cycles and amplitudes. On the other hand, 
in Paper I we showed for the first time that mi- 
crovariability is a phenomenon as common in radio 
quiet quasars (RQQs) as in core-dominated radio 
loud quasars (CRLQ), thus indicating that OM 
in quasars is independent on the radio properties, 
pointing to a discrepancy with respect to the ori- 
entation scheme. 

The orientation scheme assumes that OM is 
generated in the jet. Even if only to test this 
scheme it is relevant to investigate whether OM 
does indeed origin in the jet, the accretion disk, 
or both. Disk phenomena are mostly related with 
thermal processes such as hot spots or other disk 
instabilities, while jet emission is associated with 
non-thermal activity such as the relativistic parti- 
cle ejection events at the base of the jet. So in prin- 
ciple we can distinguish the OM origin, if we can 
identify the OM event as the result of a thermal 
or a non-thermal phenomenon. We will propose 
below a color study that can help us find some dis- 
tinctive characteristic of the variable component, 
which can be seen in the variations in the shape 
of the continuum spectrum. 

The paper is organized as follows: we begin in 
Section [2] by describing the observations, the data 
treatment, and the selection criteria. In Section [3] 
we show the results for the detected microvaria- 
tions. In Section [4] we discusses these results, and 
finally we summarize our work and give conclu- 
sions in Section [5] 

2. SAMPLE SELECTION, OBSERVA- 
TIONS AND DATA REDUCTION 

Because selection effects may yield unwanted 
biases in the study of OM, we have minimized 
them using the same selection criteria as in Pa- 
per I. Briefly, two quasar samples were observed 
and compared: one constituted by 22 RQQs, and 
a control sample conformed by 22 RLQs (see Table 
[T]). The samples were chosen in such a way that 
for each RLQ there is a RQQ with similar bright- 
ness and redshift (differing less than 10%), avoid- 
ing possible evolution biases between the sam- 
ples. To avoid orientation biases, RLQs chosen 
are all core-dominated RLQs (CRLQs), but avoid- 
ing to include objects extremely variables such 
as known blazars. When this sample was assem- 
bled, PKS 1510-089 was not considered a blazar, 



although now it is indeed. Unfortunately, it is 
hard to avoid the possibility to contaminate a 
sample of radio loud quasars by blazars, due to 
better/new observations, particularly in other fre- 
quencies. Nev ertheless, we will show in a f orth- 
commg paper faamirez et al.llin nre^Q that 
the OM detected in this object presents character- 
istics that differ of those associates to blazars. For 
this reason, we prefer to keep it in our sample. 
Choosing CRLQs instead of RLQs is a simple way 
to avoid unknown orientation effects in the sample. 
Thus, the CRLQs are expected to be as homoge- 
neous as possible from the observational point of 
view, while their intrinsic physical properties are 
a truthful representation of the RLQs population. 
All the objects in the CRLQs sample have radio 
to optical flux ratio near or above 1000, computed 
from the 408 MHz to V-ba nd, thus they belon 
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to the RLQ population after iFalcke et al 
criterion. Table [1] shows the complete sample. A 
number identifies each couple (col. 1), "a" refers 
to RQQ component and "b" to CRLQ; the object 
name (col. 2); the coordinates (J2000.0) of each 
one (cols. 3 and 4); the reported average magni- 
tude in the V band (col. 5) and the redshift (col. 
6). 

Observations were made in several seasons, be- 
tween 1998 and 2002, with four telescopes dis- 
tributed in Mexico and Spain. In Mexico, two 
telescopes, of 1.5-m and 2.1-m, were used (Mxl 
and Mx2, hereafter, respectively). These are op- 
erated by the Observatorio Astronomico National 
(OAN), located in San Pedro Martir, Baja Califor- 
nia. Observations in Spain were made using two 
telescopes. The 1-m Jacobus Kapteyn Telescope 
(JKT), located in the Roque de los Muchachos 
Observatory, and operated by the Isaac Newton 
group. We also used the 1.5-m telescope at the 
Estacion Observational dc Calar Alto (EOCA), lo- 
cated in the Astronomical Hispanic- German Cen- 
tre in Calar Alto (CAHA), which is operated by 
the Spanish Observatorio Astronomico National. 

Filters BVR of the Johnson-Cousins series were 
used throughout all the observations. The detec- 
tors used were: in Mxl, a SITE SI003, 1024x1024 
pixels of 24-/zm 2 which has a methacrome II cover 
and a VisAr, to improve the response in the blue; 
in Mx2, a Thomson TH7398 2048x2048 pixels of 
24-fim 2 which has a methacrome II cover; in JKT, 
a TEK1024AR constituted by 1024x1024 pixels 
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Tabic 1:: The sample. 







R.A. 


DEC. 






Pair 












ID 


object 


(J2000.0) 


(J2000.0) 


V 


z 


la 


US 1867 


08:53:34 


+43:49:01 


16.4 


0.513 


lb 


3C 334.0 


16:20:21 


+17:36:29 


16.4 


0.555 


2a 


1628.5+3808 


16:30:13 


+37:58:21 


16.8 


1.461 


2b 


3C 298 


14:19:07 


+06:28:35 


16.7 


1.439 


3 a 


1222+023 


12:25:17 


+02:06:56 


17 


2.05 


3b 


PKS 0421+019 


04:24:08 


+02:04:30 


17 


2.048 


la 


IE 15498+203 


15:52:02 


+20:14:02 


16.5 


0.25 


4b 


PKS 1217+02 


12:20:11 


+02:03:42 


16.5 


0.24 


5a 


US 737 


09:35:02 


+43:31:12 


16.3 


0.456 


5b 


PKS 1103-006 


11:06:31 


-00:52:52 


16.4 


0.426 


6a 


0214-033 


02:17:29 


-03:08:08.6 


16.8 


0.323 


6b 


PKS 2208-137 


22:11:24 


-13:28:09.7 


17 


0.39205 


7a 


CSO 233 


09:39:35 


+36:40:01 


17 


2.03 


7b 


PKS 1022-102 


10:24:56 


-10:31:44 


17 


2 


Sa 


CSO 21 


09:50:45 


+30:25:19 


17 


1.19 


8b 


PKS 1127-14 


11:30:07 


-14:49:27 


16.9 


1.187 


9a 


TON 156 


13:21:16 


+28:47:19 


16.6 


0.549 


9b 


PKS 1327-21 


13:30:07 


-21:42:04 


16.7 


0.528 


10a 


TON 133 


12:51:00 


+30:25:42 


17 


0.56 


10b 


3C 281 


13:07:53 


+06:42:13 


17 


0.599 


11a 


CSO 18 


09:46:36 


+32:39:51 


17 


1.3 


lib 


PKS 0514-16 


05:16:15 


-16:03:08 


16.9 


1.278 


12a 


E 0111+388 


01:13:54 


+39:07:45 


16.7 


0.234 


12b 


PKS 2247+14 


22:50:25 


+14:19:50 


16.9 


0.237 


13a 


MRK 1014 


01:59:50 


+00:23:41 


15.6 


0.163 


13b 


PKS 2349-01 


23:51:56 


-01:09:13 


15.3 


0.173 


14a 


Q0050-253 


00:52:44 


-25:06:51 


16.1 


0.626 


14b 


PKS 2243-123 


22:46:18 


-12:06:51 


16.4 


0.63 


15a 


US 3150 


02:46:51 


-00:59:31 


16.8 


0.467 


15b 


PKS 0003+15 


00:05:59 


+16:09:49 


16.4 


0.45 


16a 


US 3472 


02:59:38 


+00:37:36 


16.6 


0.532 


16b 


PKS 0122-042 


01:24:34 


-04:01:05 


17 


0.561 


17a 


MRK 830 


14:50:26 


+58:39:44.8 


16 


0.21 


17b 


OX 169 


21:43:35 


+17:43:49.5 


15.7 


0.213 


18a 


PG 1444+407 


14:46:46 


+40:35:06.0 


15.6 


0.267 


18b 


PKS 2135-14 


21:37:45 


-14:32:55.1 


15.5 


0.2 


19a 


PG 1543+489 


15:45:29 


+48:46:09 


16.4 


0.4 


19b 


PKS 1510 -08 


15:12:50 


-09:05:60.0 


16.5 


0.361 


20a 


MC3 1750+175 


17:52:46 


+17:34:20.3 


15.5 


0.507 


20b 


PKS 2128-12 


21:31:35 


-12:07:04.5 


16.1 


0.501 


21a 


PG 1538+477 


15:39:33 


+47:35:33 


16 


0.77 


21b 


PKS 1424-11 


14:27:37 


-12:03:54 


16.4 


0.805 


22a 


PG 1407+265 


14:09:23 


+26:18:18 


15.9 


0.944 


22b 


PKS 2145+06 


21:48:05 


+06:57:35 


16.4 


0.99 



Radio quiet quasars are identified with the letter a; 
radio loud quasars are identified with a letter b. 



of 24-/zm 2 was used, this chip is covered with Ar; 
and, in the EOCA telescope, a CCD Tektronics 
TK1024AB of 1024x1024 pixels of 24-/im 2 . In all 
the cases the detectors were binned to allow for a 
fast and low noise readout. 

The observational strategy was performed using 
an Analysis of Variance (ANOVA) design, similar 
to the one described in Paper I, but considering 
three filters instead of only one. Shortly, each time 
that an object was observed, it was recorded in five 
observations of ~ 1-minute exposure each, in each 
filter in the sequence BVR, except for a couple of 
cases, where the sequence was inverted to RVB. 
It was possible to monitor from 2 to 4 pairs per 
night, and 2 to 5 times each pair. Each object was 
observed at moderate air masses, always at least 
30 degrees above the horizon. The magnitude dif- 
ference errors were obtained directly from the ob- 
servations by means of the differential photome- 
try. These errors range typically from ~ 0.001 to 
0.01 instrumental magnitudes. Additionally, each 
pair was observed on the same nights and in over- 
lapping sequences, to avoid possible biases due to 
atmospheric and/or instrumental conditions. 

Previous to the extraction of data, images 
were corrected for bias and flat fields. The flat 
fields were sky flat-fields acquired at the beginning 
and/or end of every night. Objects were observed, 
insofar as possible, when they were meridian cross- 
ing, minimizing color effects by the atmospheric 
extinction. For each object, field-stars were used 
as comparison and reference stars to obtain dif- 
ferential photometric data used in the posterior 
analysis; typically, from four to eight stars were 
used (except for a couple of cases where we only 
had two stars). The selected stars were those that 
remained stable during the observations at a 20% 
confidence level. Here the results with only two 
of them will be shown. Unlike Paper I standard 
stars were observed to obtain the flux level of the 
obse rved objects as w ell. These stars are all taken 
from lLandoltl (Il992h . 

The IRAF/APPHOT0 package was used to per- 
form data reduction. The aperture radii used 
ranged from three to six arcseconds, while the 



X IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



3 



size of the sky annuli was of 10 arseconds for 
the inner radius, and 16 arcseconds for the outer 
radius. When the host galaxy was detectable, 
we took a sufficiently large aperture to avoid ef- 
fects due to the s eeing This was done although 
Carini fe Miller I (|l99ll ) and IClements fe Carini ' 



J2001I) found that the host galaxy cannot provoke 
false microvariations due to l ight that enters and 
leave s the aperture (see also iKidger fe de Diegol 
1992h . However, Cellone, Romero & Combi (2000) 
found that the host can influence spurious vari- 
ations detection, but only under extreme atmo- 
spheric conditions (with FWHM variations of ~ 5 
arcseconds) for galaxies as bright as the nucleus. 
Similar conditions were absent in our study. 

3. RESULTS 



Searching for microvariability for each object 
was done with the ANOVA, particularly, the one- 
way ANOVA test was utilized. Contrary to a 
standard test, where the object standard devi- 
ations are analyzed with respect to comparison 
stars ( e.g., a n b-i loe.nt = 2.576, for detection of var i- 
ability, Ijang fc Mille7lll997t iRomero et al~lll999h , 
ANOVA is an empirical very robust statistical 
test. We used a significance level of 0.1% (cor- 
responding to 3<r) for the one-way ANOVA test as 
variability criterion, although for some objects, in 
some filter, we will discuss also possible variations 
considering 3% and 5%. The functional procedure 
on the ANOVA test is given in the next lines. 

We have implemented One-Way ANOVA to 
search for variability comparing data measured in 
fc groups of observations sampled from a quasar 
Hghtcurve. If yij represents the value of the ith 
observation (i = 1, 2,...rij-; note that we have cho- 
sen i = 1, 2, ...5 for all the groups) on the jth 
group (j = 1, 2, ...k), we can express a mathemat- 
ical model describing each single observation: 



Vij = V + 9j + dj 
where y is the mean of the whole dataset, gj is the 
between-groups deviation (gj = y~j — y) and 
the within-groups deviation, often called residual 
or measurement error (e^ = yij — yi). 

ANOVA tests whether the means of the groups 
are equal. Statistically, this can be expressed say- 
ing that the null hypothesis (i.e., the hypothesis 
that we are testing) is that the means of the dif- 
ferent groups are the same, and the alternate hy- 



pothesis (i.e., the hypothesis that we will accept if 
the null hypothesis would be rejected) is that at 
least one group has a mean different from the oth- 
ers. In our case, the alternate hypothesis implies 
a microvariability detection. 

Following the implicit mathematical model of 
ANOVA, the total sample variation can be split 
into variation between and within groups: 

k rij k k rij 

j=l i=l j = l j=l i=l 

where the left term is a measurement of the total 
deviations of the data with respect to the mean, 
the second term measures the total variation be- 
tween groups, and the third term the total errors. 
The equation can be abbreviated to: 



55 



Total 



S Si 



Group 



SSf 



When the null hypothesis is true, the k groups 
sample data will be normally and independently 
distributed, with mean /i and variance a 2 . Thus 
the statistic: 



F = 



MS G 



roup 



MSt 



S S Group / (k 1) 
SS E rror/(N - fc) ' 



will follow an F distribution with k — 1 and 
n — k degrees of freedom, and where the pseudo 
variances MScroup and MSError are mean esti- 
mates for the variations between groups and er- 
rors, respectively. Given a certain significance 
level a, if the F statistic exceeds the critical value 
F(k-x,n-k,ot) t ne null hypothesis should be re- 
jected. 

The rest of this section is devoted to describe 
the results, object by object, for those cases where 
variability was detected. In Figures Q] we show the 
light curves for differential photometry for each 
object. Symbols represent the mean for each group 
of five observations of each object (left panel) and 
the stars (right panel). Error bars correspond to 
a standard error toward each side, and they were 
calculated from the dispersion of each group. 

From a total of 130 observations, we detected 
9 microvariability events, and 2 possible ones, in 
eight out of 44 observed objects, three CRLQs and 
five RQQs. Two quasars displayed microvariabil- 
ity in more than one night: the CRLQ 3C 281 was 
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observed to vary in two nights, while microvari- 
ability in Mrk 830 was detected one night, and 
probably were present in two more occasions. 

3.1. Radio Loud Quasars 

3C 281.- This radio loud quasar is known to 



resid e in a very rich galaxy cluster (jYee fc Green 



1987t h A variation of 0.15 magnitude s in one year 



has b een reported in the infrared (jEnva et al 



2002al lbf). We observed this object for Paper I 



on May of 1997, but microvariability was not de- 
tected. We observed this quasar with EOCA on 
2000 March 1 and 4, detecting OM events in both 
nights. 

2000 March 1.- The object was observed dur- 
ing 2.7 hours three times. A visual analysis of the 
light curves shows a decrease in brightness of al- 
most 4 hundredths of a magnitude (Figure 1(a)). 
AN OVA confirms a microvariation in B and V, de- 
tectable at a 0.1% significance level, but not in R. 
This detection is appreciated with all the available 
stars. In V, the amplitude is AV = 0.031 ± 0.007 
mags, in 2.7 hrs, while between the second and 
third groups AV = 0.024 ± 0.007 mags, in -1 
hour. In B, the amplitudes arc AB = 0.05 ± 0.01 
mags between the first and the third observations 
and AB = 0.03 ± 0.01 mags between the second 
and the third. 

2000 March 4.- 3C 281 was observed during 
3.75 hours in three occasions. A visual revision 
of the light-curves establishes the occurrence of 
a possible variation in the three bands (Figure 
1(b)). However, ANOVA only detected a varia- 



tion in V and a marginal one in R (with a sig- 
nificance level of 5%). The variation amplitude is 
AV = 0.027 ± 0.006 mags between the first group 
and the third. 

PKS 1510-089.- This RLQ has shown in- 
tense ac tivity on long and inte r mediate-term vari 



abilit y (iLiller fc Lillerl Il975t iTornikoski et al 



11994 IVillata et al. I Il997l ), a s well as corre 



lated optical/radio variability (ITornikoski et al 



1994 ) . OM also has been reported (|Villata et al 



1997; IXie et~aT~l l200ll;lDai et al 



20041 ). In particular. iDai et al 



2001; XieetaL 



200 if ) reported a 



variation of 2 mag in 41 mi nutes in R band at the 



(1999) 



end of 2000 May. Howeve r, iRomero et al. 
and Romero et al. I ( 2002 ) did not detect any ac- 



tivity when they observed it between 1998 and 



> 0,94- 



I i 
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(a) 3C 281, 2000 March 1 



02:00 03:00 04:00 05:00 
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02:00 03:00 04:00 05:00 

U.T. 



(b) 3C 281, 2000 March 4 

Fig. 1. — : Light curves for the objects discussed 
in the text. In each Figure, left panel shows the 
data for the quasar, while right panel for the com- 
parison stars. Symbols represent average values of 
five consecutive observations, and the error bars 
have been estimated from data dispersion. Foot- 
note of each Figure indicates the quasar name and 
the observation date. 



1999. Our observations show that PKS 1510-089 
was in a stage of activity during the first semester 
of 2000. We monitored this object during 4.3 
hours in 5 different occasions during the night of 
2000 M arch 20 (two mon ths before the observa- 



tions of iDai et al. 1 120011 ). Light-curves (Figure 



1 (c) I indicate the presence of an OM event de- 
tected in V and R. While a progressive increase in 
brightness is detected in V, a variation is observed 
between the first and the second datasets in R, 
which remains in the same level of brightness the 
rest of the session. ANOVA agrees with the visual 



revision of Figure 1(c) This variation is confirmed 



with a 0.1% level of significance. The comparison 
between the stars indicates that they do not vary. 
In V, the microvariation seems to evolve during 
the 4.3 h of observation, with a total amplitude 
of 0.104 ± 0.014 mags (with respect to the first 
data set changes are: AV^ = 0.03 ± 0.02 mags in 
1 hour, AV = 0.07 ± 0.01 mags in 2.2 hours, and 
Ay = 0.09 ± 0.01 mags in 3.25 hours). In R, the 
amplitude between the first dataset and the aver- 
age of the subsequents is AR — 0.108±0.008 mags. 
From the light-curves, it is clear that color changes 
occurred during this OM (see Figure 1(c)). 

PKS 0003+15.- Historically, this quasar 
has shown intense activity on long-term and 
intermediate-term variability (iBarbieri et al. 119791: 



1995, 1998: Garcia et al 



Pica et al. ll988l:ISchramm et al Hl994l : lGuibin et al 



19991 ). 



OM has also 



been detected. On October of 1994, Jang fc Miller 



([19951 ) observed a variation in R with an ampli- 



tude o f 0.11 mags in 2 hours, while lEggers et al 



(200 0J) report e d cha nges of 0.005 mag/h in R. 
Garcia et all (|l999h also reported microvariabil- 
ity for this quasar. PKS 0003+15 was reported 
in Paper I with no evidence of microvariability. 
We have monitored this quasar during 3.13 hours 
in four occasions on October 11 2001. Figure 



1(d) shows that while brightness remained con- 
stant in B and R, a decrement in V was de- 
tected. ANOVA confirms this visual revision. The 
brightness for this object showed a decrease of 
AV^ = 0.025 ± 0.005 mags during the monitoring 
time. 

3.2. Radio Quiet Quasars 

MC3 1750+175.- The only variability event 
previously reporte d for this quasar corresponds to 
radio wavelengths ( Fanti et al. Ill98ll iMantovani 
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(d) PKS 0003+15 
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19821 ). On 1999 August 20, we observed this ob- 
ject three times during 3.13 hours. The light- 
curves show a microvariability event only visible 
in the V band (Figure |l(e)[ ). ANOVA confirms 
this appreciation, the amplitude of this variation is 
AV" = 0.04 ± 0.008 mags between the first dataset 
and the third. 

CSO 21.- Observations for CSO 21 were re- 
ported in Paper I, but without a positive detec- 
tion of microvariability. For the present study, the 
object was observed in three occasions on 2001 
December, during 2.5 hours. A gradual brightness 
increase is appreciated in the three bands (Fig- 
ure 1(f) I. However, ANOVA indicates a varia- 
tion in V only, a barely marginal variation in R 
(with a significance level of 10%), and no varia- 
tion in B (if the observed tendency in the B light- 
curve is a variation, it is hidden by the internal 
data dispersion). The amplitude of the variation 
is AV = 0.025±0.005 mags in 2.5 hours. Variation 
amplitude of the R band is AR = 0.016 ± 0.004 
mags between the first dataset and the third. 

1628.5+3808.- In Paper I we reported the 
cases of microvariability for this object for the first 
time. They are two events on 1997 February an 
May. In February, a variation with amplitude of 
0.082 mags throughout 2.4 hours; while in May, 
throughout six monitoring hours, the OM con- 
sisted of pulsations with amplitudes from 0.04 to 
0.057 mags. We observed this quasar in two occa- 
sions on 2001 June 13. The light curves show an 
appreciable brightness decrease in the three bands 
(Figure 1(g)). ANOVA indicates that this varia- 



tion is detectable in the B and R, but in V the 
variation is only marginal (the significance level is 
3%). The data dispersion in the comparison stars 
do not show evidence for statistical differences. 
The OM amplitudes are AB = 0.061 ±0.011 mags 
and AR = 0.018 ± 0.004 mags, in the 2.39 hours 
of monitoring. The V marginal detection would 
have an amplitude of 0.03 ± 0.01 mags. 

US 3472.- In Paper I we reported an OM event 
for this quasar with an increase of 0.02 mags in 50 
minutes. During 2 hours, a series of oscillations 
followed to this brightness increase, with ampli- 
tudes between 0.040 and 0.055 mags. Long-term 
optical and infrared variab ility has also been re- 
ported ( Enva et al. Il2002ah . We observed this ob- 
ject during 2.5 hours in three occasions on 2001 
December 20. The light-curves show a brightness 
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decrease in l^and R (Figure 1(h)) 
firms the variation 



ANOVA con- 



In V, the microvariation seems 
to have kept a constant rate with a total ampli- 
tude of AV = 0.020 ± 0.003 mags in 2.35 hours. 
Between the first and the third dataset in the R 
band, the brightness varied by 0.014 ±0.002 mags. 

Mrk 830.- Long- and short-term variability, 
and microvaribility have been reported for this 
quasar. We observed this object in four differ- 
ent occasions, OM events being detected on 1999 
August 17 and 18, and on 2001 June 14. The 
object possessed a brightness very similar on the 
three occasions, V ~ 17. 6. This brightness issim- 



ilar to that reported b y IStenanian et al. I (2001) 
(V = 17.29±0.03), and lChavushvan et al. I (jl995l ) 
(V = 17.62). Unfortunately, the fields of the two 
telescopes used are different, and thus it was more 
convenient to use a different set of comparison 
stars for the August of 1999 and the June of 2001 
observations. 

1999 August 17.- The object was observed dur- 
ing 2.13 hours in three times. A visual revision of 
the light-curves indicates the possible detection of 
a microvariability event in V (Figure l(i)l. Due 
to instrumental problems, in the first sequence 
of the R band, we missed the object for several 
images. Thus we can present only the complete 
curve for the comparison stars. ANOVA shows 
that this variation would be marginal, with a sig- 
nificance level of 3%. The increment in brightness, 
if real, has a constant rate with a total amplitude 
of AV = 0.03 ± 0.01 mags in 2.31 hours. 

1999 August 18.- The monitoring dataset con- 
sisted on three observation groups taken during 
two hours. Again, a revision of the light-curves 



shows a possible variation in V (Figure l(j)l. 



ANOVA detects this feature as a marginal vari- 
ation, with a significance level of 3%. As on the 
previous night, the tendency is a variation at a 
constant rate, but in this time with a total ampli- 
tude of Ay = 0.05 ± 0.03 mags. 

2001 June 14-- On this night the quasar was 
observed twice. A visual revision of the light- 
curves indicates a possible variation in the three 
bands (Figure 1 (k) ) . The brightness increase is 
a few hundredth of magnitude between the two 
groups of observations. ANOVA confirms the re- 
liability of this variation. While in V and R the 
significance level is 0.1%, in B the value is ~ 5%. 
The amplitudes are AV = 0.039 ±0.007 mags and 



8 



-0,3'rf 
-0,38 
-0,37 

E 

1 -0,36 
^ -0.35 
-0,34 

-0,46 
-0,45 
E -0,44 
.E -0,43 

> 

< 

-0,42 
-0,41 
-0,49 
-0,48 
E -0,47 
| -0,46 
< -0,45 
-0,44 



21:00 22:00 23:00 

U.T. 



21:00 22:00 23:00 

U.T. 



(i) Mrk 830, August 17 1999 



-0.42 
-0,40 

-0,38 

E 

~ -0,36 
< -0,34 
-0,32 

-0,48 
-0,46 
p -0,44 

>~ -0,42 
-0,40 
-0,38 
-0,50 
-0,48 
2 -0,46 

£ -0,44 
-0,42 



1 I 



I I 



22 00 

U.T. 



22:00 

U.T. 



(j) Mrk 830, August 18 1999 



Fig. 1. 



Continued. 



3,12 

3,14 

E 

1 3 ' 16 
CD" 

< 3,18 



3,08 
3,10 



C£ 3,12 
< 

3,14 



1,14 

| 1,16 

DCf 1,18 
< 

1,20 

1,48 
1,50 

| \» 

5 LSI 

1,56 

1,70 

1,72 

I 1.74 

OZ" 1,76 
< 

1,78 



07:00 07:30 

U.T. 



07:00 07:30 08:00 

U.T. 



(k) Mrk 830, June 14 2001 



Fig. 1. 



Continued. 



AR = 0.030 ± 0.007 mags, occurring during 1.18 
hours. If a variation exists in B, the amplitude 
would be of AB ~ 0.04 mags. 

4. DISCUSSION 

An immediate result of this work is that mi- 
crovariability is a more complex phenomenon than 
appreciable in single optical band studies. The 
different color evolutions during microvariability 
events will be analyzed in a forthcoming paper 
( Ramirez et al.llin preparation ) . In this section we 
will compare statistically the results of OM ob- 
tained from the two samples of RQQs and CRLQs. 
The aims are twofold: 1) to find inferences on the 
populations of which the samples have been ex- 
tracted when comparing each pair of objects, and 
2) to describe the probability of detecting a mi- 
crovariability event as a function of the duration 
of the observation for each object. Several well- 
known statistical tests are applied to attain these 
objectives. Some of these tests deal with the qual- 
itative properties of the results, such as detection 
or non-detection of OM (§ 14. ip . and others deal 
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directly with the numerical results (§ I4.2[) . Be- 
cause of the differences in the variability detected 
for each band, the numeric tests were carried out 
separately. 

Statistical tests are described in detail in Pa- 
per I. Shortly, a \ 2 test for homogeneity is used 
to compare the occurrence of OM events between 
the CRLQs and the RQQs samples. Using this 
test, one can determine if a quasar type presents a 
statistically larger number of OM events than the 
other. On the other hand, a quantitative test on 
the differences among the samples is carried out 
through the comparison of the variances of the 
paired objects. The observed variance was con- 
sidered to arise from two different error sources: 
one corresponding to the observational error and 
another due to the intrinsic variance of each ob- 
ject (see Paper I for details). Then a couple of 
statistical tests were applied to the mean of the 
differences of intrinsic variances of RQQs and CR- 
LQs: a t- Student matched pair test and an inde- 
pendent sample test. At the end of the section, the 
microvariability occurrence as a function of moni- 
toring time, i.e., the duty cycle will be analyzed. 



The variations showed in Figs. 1(c) and 1(d) 
for PKS 1510-089 and PKS 0003+15, respectively, 
differ very much in different bands. However, 
these variations can be explained in the con- 
text of emission models. We wil l explain such 



variations in a forthcoming paper (jRamirez et al 
in preparation ). 



4.1. x 2 Test for Homogeneity 

The \ 2 test for homogeneity is used to compare 
two or more qualitative properties of several sam- 
ples when the number of elements in each sample 
is fixed. In our case, we will split the observations 
of the two samples of RQQs and CRLQs into those 
that show variability and those that do not. 

For each sample of RQ and CRL objects the 
microvariability events were counted regardless of 
the band of detection. As in Paper I we note 
that although each observation will be considered 
independently, some of them may belong to the 
same object. In doing this, some unwanted effects 
can be introduced because the observations of the 
same object may not be completely independent. 
This may be the case, for example, if one object 
is intrinsically much more variable than another. 



The total number of paired observations for 
each sample is 63 (although the number differs for 
each band, i.e., the occasions in which one could 
observe the two members of the same pair were 61 
for B, 63 for Fand 61 for R). In 9 of these monitor- 
ing microvariability was detected for at least one 
member of each pair, with a significance level of 
0.1%. Five out of these 9 were detected in RQQs 
and four in RLQs. This number rises to 11 when 
we take into account the two marginal detections 
of Mrk 830. 

Table [2] shows the x 2 contingencies for these 
data. Column 1 displays the number of OM events 
with a significance level of up to 0.1%. The num- 
ber of observations in which there are no evi- 
dences of variability appear in column 2, the total 
number of observations appear in column 3, and 
the expected frequencies of microvariations under 
the null hypothesis, i.e., if both types of objects 
possess the same microvariability properties, are 
shown in columns 4 and 5. These last columns are 
calculated multiplying the partial totals of each 
line and each column and then dividing by the to- 
tal number of observations. The homogeneity test 
shows that the probability of observing the num- 
bers in Table [5] can be obtained 27% of the times 
when extracting two samples in an aleatory way 
from the same population. The probability to ob- 
tain the observed OM frequencies is large enough 
and thus the null hypothesis cannot be rejected. 
In other words, it is reliable to assume similarity 
between the two samples. 

Table 2:: x 2 contingencies table. 





Observed 




Prospective 


Quasars 


V 


NV 


Total 


V NV 


RQQ 


5 


58 


63 


4.5 58.5 


CRLQ 


•1 


59 


63 


4.5 58.5 


Total 


9 


117 


126 


9 117 



Contingencies Table to compare the radio properties of the 
microvariability of the two quasars samples. It does not 
seem to have any difference when the three niters are con- 
sidered. This result indicates that optical microvariability 
is a common phenomenon to the object classes, a = 0.01 
as variability criterion. 



10 



4.2. Quantitative Test for Differences be- 
tween RQQs and CRLQs 

Although the test for homogeneity did not show 
differences for the occurrences of OM between 
RQQs and CRLQs. a test on the quantitative val- 
ues of the variance for both samples of objects 
might show whether the OM is larger in one type 
than in the other. The color changes during vari- 
ations make it interesting to consider observations 
of each band separately. 

In Paper I we showed that it is better to es- 
timate the difference between RQQs and CRLQs 
breaking the variance into two components than 
considering a single observed variance. The in- 
trinsic variance for a source, V r , can be written as 
V r = V — e 2 /5, where V refers to the observed 
variance and e 2 are the squared errors of the data 
weighted by the number of observations in each 
group (five). The last refers to the inherent error 
of the observation and reduction data processes, 
and it is calculated empirically from the variance 
for each of the five observation groups. 

The t-Student test can be used to find differ- 
ences between the means when these have a nor- 
mal distribution. This is approximately the case 
for the intrinsic variances difference. For these 
data, we have computed the matched pair test, in 
which the mean of the distribution of these dif- 
ference should not be significantly different from 
zero if both samples are extracted from the same 
population. The test indicates that there are no 
significant differences among the samples: it gives 
a probability that the observed differences occur 
by chance of 18.7% for B, 19.1% for V and, 87% 
for R. When an independent pair test is carried 
out, the mean difference values are similar, within 
the errors, and so are the uncertainties of such dif- 
ferences. 

This test was repeated considering average val- 
ues of the intrinsic variances for each pair of ob- 
jects. The qualitative result remains the same, 
however quantitative results do change: in V, the 
probability that the observed differences happen 
by chance increases to 28%, in B changes to 23% 
and in R it diminishes to 53%. The averages in this 
case are distributed so that the fits are quite good, 
because although we have lost degrees of freedom, 
the data of the variations themselves are masked 
when they are averaged with other data. It is nec- 
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Fig. 2. — : Histogram for the intrinsic variance 
differences among CRLQ and RQQ paired. The 
continuous line represents the result of adjust- 
ing the mixture of two normal distributions to 
the data without detected microvariability (dot- 
ted line) and those with detected microvariability 
at least in a single band (dashed line). 



essary to note that the probability obtained for the 
data of the Fband is consistent with the value ob- 
tained in Paper I, where the same filter was used. 
These results imply that the incidence of OM is 
independent from luminosity and redshift. 

The tests neatly point out to an extremely small 
and non significant difference in the microvariabil- 
ity displayed by the two samples, yielding no evi- 
dence that they come from different populations. 
Figures 2(a)[ |2(b)| and 2(c) show the histograms 
for the intrinsic variance differences of the obser- 
vations paired for each filter. Except in the case of 
R, for which data can be fitted using a single Gaus- 
sian profile with standard deviation of 2.3 x 10~ 3 
and an average value of ~ 0, the others cannot 
be accurately fitted by a single Gaussian. Taking 
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Fig. 2. — : Continued. 



into account the results presented in Paper I, we 
separated the data into observations with (dashed 
line) and without (dotted line) detected variations. 
Then, good fits are obtained. As in Paper I we 
found that some variations could occur and not 
been detected. An OM event could occur in the 
three bands even if it is detected in only one. 

In the B band, the data group where OM was 
not detected is centered in ~ 0, and it has a stan- 
dard deviation of 4 x 10 -3 , i.e., indistinguishable 



from zero (Figure 2(a) ) while, the group with vari- 
ability has a standard deviation of 4 x 10~ 3 and is 
centered at 1 x 10 -3 . These results indicate that 
the samples are not statistically distinguishable. 



In V there is a similar behavior (Figure 2(b) I. The 



non variable component was fitted with a Gaus- 
sian that has a standard deviation of 1.2 x 10~ 3 
and is centered on 4 x 10 -4 , while the variable 
component is centered at ~ and has a standard 
deviation of 4 x 10~ 3 . The data with OM detec- 
tion were fitted by a Gaussian centered on 4 x 10 
with standard deviation 3 x 10 -3 while the non 



detection data were best fitted with a Gaussian 
centered on with standard deviation 2.3 x 10 -3 
(Figure 2(c)). These results can change when we 



consider that some quasars were observed in more 
occasions than others. 

All the statistics discussed here indicate that 
there are no significant differences in the OM prop- 
erties between CRLQ and RQQ. While this is at 
odds with some other studies, it is consistent with 



the re sults rep orted in Paper I by IStalin et al. 
(<2004h . and bv lGupta fc Joshi I (|2005t) . 



4.3. Sample Selection and Duty Cycle 

We can characterize the probability with which 
an OM event is detected as a function of the dura- 
tion of the monitoring. The Duty Cycle (DC) is a 
very useful tool to compare the behavior of differ- 
ent types of AGN, and to plan future observations. 
The DC is defined as the ratio of the number of 
observations with a positive detection of OM to 
the total number of observations, weighted by the 
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period of tim e during which the ob jects have been 



observed (see Romero et al. ~l ll999l) . Of course OM 



depends also on intrinsic causes that determine its 
occurrence. The last however, are unpredictable 
and unquantifiable. OM is believed to be a tran- 
sient phenomenon. Independently of its origin, 
neither thermal nor non-thermal instabilities arc 
expected to be permanent. We have calculated 
the D C using expression (2) from iRomero et al. I 
( 19991 ). but considering the joined observations for 
RQQs and CRLQs, because we detect no differ- 
ence between both samples. When the DC is cal- 
culated without taking into account the band, a 
value of 8.45% is obtained. If marginal detections 
are considered, the value does not changes too 
much: 10%. Separating data by band, we obtain 
DC B = 2.3%; DC V = 6.8%, and DC R = 5.8%. 

The DCy is three times lower than the one re- 
ported in Paper I, where the same filter was used. 
The discrepancy could be due to the use of three 
filters in the present study. When we analyze the 
light curves from Paper I, we can see that one 
third of the events correspond to variations with 
constant increase or decrease. The rest correspond 
to changes with oscillating behavior. If this pro- 
portion is a general behavior, and since we have 
now a lower temporal covering, we expect to de- 
tect only the constant changes, i.e., a third of all 
the possible detections with one filter. 



(2004) 



re- 



For the blazar class, iSagar et al. 
ported a DC = 72%, and iGupta fc Joshi 
reporte d DC = 60%, for microvariability; while, 



( 200511 



Carini I (Il990h reported a DC of - 80% for 
timescales from ho urs to days. For the quasars 
Stalin et all (|2004h calculated a DC RQQ = 17%, 
and DCrlq = 15%. Their observations were car- 
ried out in the R band. Using the data from table 



4 of P aper I, and the expression of IRomero et al 
(|l999h . a DC of 18% is obtained for the join of 
RQQs and CRLQs samples. 

Other groups have fo und differences i n the Duty 
Cycle of these objects. IRomero et al. ( 19991 ) cal- 
culated a DC of 68% and 6.9% for RL and RQ 
objects, respectively (their observations were car- 
ried out with a V fi lter). Taking the data of 
Jang fc Miller! (|l997h . the DC^qq = 7.6% (their 



observations were carried out with a R filter) 
while DCrlq = 60.6%. With the data from 



ICopal-Krishna et al. I (l2000l ). DC 



J RQQ 



6.1% 



least are very probable, while when including those 
probable detections T)Grqq = 29.4%). These re- 
sults show a clear difference between radio loud 
and radio quiet quasars, which represents a behav- 
ior very different to the one found in the present 
study. 

At this point we want to stress the relevance 
of a good sample selection. However, all these 
groups included BL Lac objects among the RLQ 
samples, and it is very possible that this can al- 
ter the value of the DC for the RL sample. The 
DC for the RQQs on the other hand are very 
similar (~ 6 — 7%) to that found by us (for the 
complete samp le). To illustr a te th is point, we 



point out that IRomero et al. ( 1999f ) 



use a sam- 



ple of RL objects that includes several BL Lac ob- 
jects. 53% of their observations and 67% of their 
microvariability detections, correspond to blazar 
type objects. On the other hand, the RQQs sam- 
ple of Jang fc Miller 1 (1997) is constituted mainly 
by non OM objects, since it is conformed by 90% 
of low brightness quasars. The average magni- 
tudes for RQs are similar to the average magni- 
tudes for RLs (averages only differ in 2%), but 
the average redshifts are much smaller for RQs 
(average redshifts differ in 62%; for RQQ z aver- 
age is 0.5 while z average for RLQ is 1.3). Al- 
though in Paper I we have not found evidence 
of dependence of OM on brightness and/or red- 
shift, it would be necessary to investigate further 
such possible biases. The non homogeneous sam- 
ples problem could also be p resent in the work of 



Gopal-Krishna et al. 



(2000). These authors re- 
port observations of RQQ, without a control sam- 
ple, comparing with their previous results on radio 
loud objects. 

It is also worth noticing that, when the above 
authors make a sub-selection of objects with the 
selection criteria described here and in Paper I (ex- 
cept for the use of lobe dominated radio quasars 
instead of CRLQ), they obtain similar compara- 
tive results to thos e reported in Paper I (see also 
Stalin et al. Il2004l) . 



(not taking into account their detections that at 



Finally, in Paper I we developed an alterna- 
tive method to calculate the Duty Cycle. This 
method considers that all the quasars present the 
same DC in H hours of monitoring. We grouped 
the data in parcels according to monitoring du- 
ration, with intervals of 1 hour. Parcels larger 
than 6 hours have been excluded from the sam- 
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pling since they were very few for a good statistics. 
Using the current data in the equations 1, 2, and 
3 from Paper I, we get that the mean value for the 
probability of detecting an OM event in 1 hour of 
monitoring, Pi has a value of 8.9 ± 2%, which is 
marginally different from the value found in Paper 
I (Pi = 5 ± 2%, however take also account the fact 
commented above that in this study the detection 
of oscillating variability may be dumped by the 
larger lags in the monitoring). Table [3] is similar 
to Table 4 in Paper I. Columns show, respectively, 
the values of the number of monitoring hours, H; 
the number of monitors in H hours, Nh', the num- 
ber of OM events detected in H hours, Eh] the 
probability of detect an OM event in H monitor- 
ing hours, Pjj'i and the probability of detection in 
1 monitoring hour, P\. 

4.4. ANOVA and X 2 Tests 

To illustrate the reliability in the use of 
ANOVA, additionally we have applied a test x 2 
to the data. The average, P>i, has been taken 
from each group of five observations, and the er- 
ror, Si, is obtained from the dispersion of these 
five observations. Table [4] shows the differential 
photometric data for each group of observations 
(ti,t 2 and t 3 of US 3472). Taking the whole data 
set mean, < B >, we can obtain the statistical 

2 _ x ^ (Bj- < B >) 2 

X s 2 

As example, we see the case of US 3472 for 
the V band. Please note that we are carrying a 
large number of decimal places during the calcula- 
tions. For this quasar we have the following values 
obtained from differential photometry, comparing 
with a stars field 

Bi = 0.5602 ± 0.00206, 



B 2 = 0.5718 ±0.00116, 
B 3 = 0.5808 ±0.0022. 

From these values, the median is < B >= 
0.57093. Thus, we have x 2 = 27.13095 + 0.5625 + 
20.12746 = 47.82091. As the tabulated value for 
X 2 , with 2 degrees of freedom, for a level of mean- 
ing to the 0.1% has a critical value of 13.8150, we 
have detected variations at this level. 

For the case of stars field, the raw differential 
photometric data is shown in Table [5l and the 
means are: 

B l = -0.9446 ±0.00558, 

B 2 = -0.9422 ±0.00256, 

B 3 = -0.9466 ±0.00216, 

which median is < B >= -0.94447. With this, 
we have X 2 = 0.000542773 + 0.78627 + 0.97242 = 
1.75923. 

For comparison, using ANOVA the sum of the 
differences of the averages of the data shown in 
Tabled] is 0.00107, while the square difference of 
the internal dispersion is 0.0002084. As the de- 
grees of freedom are k — 1 = 2 and TV — k = 12, 
we have F = 30.80614203. With a level of signifi- 
cance of 0.001, a value of 12.97 is obtained. With 
such result, we have a positive detection of OM. 

With the same procedure for the field stars (Ta- 
ble HJ), we obtain that the sum of the square dif- 
ference of the averages is 0.0000485335, while the 
difference of the internal dispersion is 0.0008472. 
Thus, we have that F = 0.34721671. 

This example illustrates that microvariability 
detections with ANOVA are at least as reliable 
as those that can be obtained using x 2 ■ I n our 
example, we could estimate errors for each ob- 
servation in the x 2 ^ es ^ from the same quasar 
dataset, because the observations were performed 



Table 3:: Duty Cycle 



Table 4:: US 3472 minus star 3 data 



H 


N H 


E H 


Ph 


Pi 


h 


*2 


*3 


1 


9 


1 


0.111 


0.111 


0.558 


0.571 


0.583 


2 


21 


1 


0.048 


0.024 


0.564 


0.570 


0.585 


3 


53 


5 


0.094 


0.032 


0.555 


0.569 


0.585 


4 


25 


2 


0.080 


0.021 


0.558 


0.575 


0.576 


5 


17 









0.566 


0.574 


0.575 
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using an ANOVA design. This is not usually the 
case with error estimates in \ 2 based designs. In 
X 2 tests, errors are usually calculated either from 
differ ent data sets (one or more comparison stars) 
(e.g., Romero et al. |[T999() or from IRAF theoret- 
ical errors multipli ed by a correction factor (e.g., 



IStalin et al. 2005 ). Besides, this estimated er- 
ror is often considered as a constant along the 
dataset, missing the internal variability of the data 
due to short timescales factors (tiny clouds, seeing 
changes, etc.). By contrast, ANOVA uses inter- 
nally consistent errors calculated from the same 
dataset. 

5. SUMMARY AND CONCLUSIONS 

We present a comparative study of the occur- 
rence of OM in radio loud and radio quiet quasars 
using BVR bands. Following Paper I we stress 
the importance of sample selection criteria and ob- 
servational methodology. Microvariability events 
were analyzed using one-way ANOVA test, rein- 
forcing with a visual revision of the light-curves, 
and a \ 2 test. Five events were found in RQQs and 
four in CRLQs, with a significance level of 0.1% 
for a positive detection. Then, microvariability 
properties for each sample were compared using a 
X 2 test for homogeneity. No significant difference 
among the samples was found, in agreement with 
the results obtained in Paper I. 

Additionally, it was shown that both quasar 
types present the same quantitative microvariabil- 
ity properties. This is evidenced applying a t- 
Student test to the variance differences of each 
CRLQ-RQQ paired couple. Since quasars were 
paired by brightness and redshift, these proper- 
ties do not influence the result. Moreover, these 
results do not differ when applying a statistical t- 
Student test for independence on the total obser- 
vations, confirming the results reported in Paper 
I, i.e., that the evolution of the objects seems to 

Table 5:: US 3472: star 3 minus star 6 data 



tl 


*2 


h 


0.952 


0.944 


0.948 


0.959 


0.933 


0.939 


0.936 


0.948 


0.951 


0.928 


0.945 


0.95 


0.948 


0.941 


0.945 



have little influence on microvariability. 

Similar results have been found f or long-term 



and intermediate-term var iability (e.g. lWebb fc Malkan 



20001: IVanden et aT~ll2004f) . The union of all these 



results, and the discovery of radio jets in RQQs 



20011; iBlundell et al. 1 120031 ) suggest that the OM 



(IBlundell fe Beaslev I [1991 IBlundell fc Rawlings 



source must be the same in both types of quasars 
or at least very similar. This suggestion is re- 
inforced by the detection of variations in ra 



dio frequencies for RQQ (e.g..lFanti et al. 
Mantovani Ifl982l ; IBlundell et al. 1120031 ). 



1981 



Finally the Duty Cycle (DC; the percentage of 
monitoring time that an AGN shows variability) 
was calculated. This probability is 2.3% in the 
B band, 6.8% in the V band, and 5.8% in the R 
band. Without separating data by filter, DC= 
8.45% is obtained. These differences of DC be- 
tween bands might arise from spectral variability 
during the MO events. 

Res ults obtained by Ijang fc Miller I dl995l 



997h iGopal-Krishna et al. I (|2000h . and lRomero et al. 
( 19991 ) indicate that microvariability could be re- 
lated to radio emission properties of the quasars. 
However, these results are probably biased due 
to sample selection criteria (see 4.3). On the 
other hand, it has been shown that with selec- 
ti on criteria similar t o th ose discussed in Paper 
I. IStalin et al~l (|2004h and lGupta fc Joshi I (|2005l ) 
obtain compatible results with those of the present 
paper. 

As our results seem to indicate that the opti- 
cal microvariability does not depend on the radio- 
properties of quasars, we would not exclude the 
possibility that it can have a similar origin in both 
RLQs and RQQs (at least for the commonly de- 
fined ranges). Although this would be an impor- 
tant result, deriving from comparative studies, we 
would still have the additional problem to estab- 
lish unambiguously the origin of the microvariabil- 
ity. The possibil ity of a symbiosis disk-jet phe- 
nomenology (e.g.. lFalcke et al~ 19961 ) can compli- 
cate any analysis, since instabilities (that produce 
variability) originated in t he inner disk can be 
propagated to the jet (e.g., Gupta fc Joshi I 20051 : 
IWiitall2005l) . 

In a forthcoming paper ( Ramirez et aljin preparation ) 
we shall address this issue using our multiband 
data. Events arising in the accretion disk are 
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expected to show a thermal signature, while jet 
perturbation events should be dominated by non- 
thermal processes. Thus the spectral evolution 
during a change of brightness, in other words, a 
color variability analysis can allow us to differen- 
tiate thermal from non-thermal origins. 
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